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ABSTRACT 



The time-optimal control functions for sample control systems were 
computed using the iterative computational procedure proposed by L. W. 
Neustadt as adapted for differential analyzer solution. The computational 
procedure investigated was developed from the method of steepest ascent. 
Adaptation for differential analyzer solution required that finite 
sized steps be taken during the ascent vice the infinitely small sized 
steps permissible in theory. This presented one of the problems of the 
computational procedure. No optimization of the step size was attempted 
in this work, however, several criteria for the selection of the step 
size were used with success on the two-dimensional systems. 

The optimal control functions for two-dimensional systems were 
readily computed using Neustadt' s iterative computational procedure. 
Analytic work was performed to verify some computer solutions. Agreement 
of computer and analytic solutions was favorable. The optimal control 
function for the three-dimensional system investigated was not obtained. 
After sixty iterations th re was no apparent convergence to an optimal 
control function. Some refinement of the computational method used in 
this work would be required to extend the problem solutions to higher 
order systems. No work was performed to determine the effect computer 
errors would have on problem solutions. 
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I . INTRODUCTION 



There has been an abundance of literature written concerning the 
classical control problem, i. e., given a control system, assumed to be 
normal, such that 
(1-1) x= A(t) x + B(t) u 

where: x represents an n-dimension state vector 

A(t) represents an n x n matrix, not necessarily constant 
B(t) represents an n x r matrix, not necessarily constant, and 
u represents the r-dimension control vector 
what is the optimal control, u°, which will drive any arbitrary state of 
the system, x^tg), to some specified final state, x^t^), in the minimum 
time. It is well known that the desired control law is bang-bang and 
is of the form 

(1-2) u = sgn ( B T ( t ) X T (t 0 ,t) \?) 

where: yj represents an n-dimension constant vector 

X(t ,t) represents an n x n matrix satisfying: 

X = A(t) X , where X(tQ,tQ) = I 

Assuming || || ^ 1, where u^ are the components of u, then let 

oi = B T (t) X T (t 0 ,t)^ , then (sgn 5* ) j _ = +1, o< > 0, (sgn cK ) i = -1, 
o< i < 0, and (sgn c< is undefined for * 0. 

Lucien Neustadt published a method for the synthesis of the control 
function which was adaptable to analog computer and digital computer solution, 
refer to L. W. Neustadt, "Synthesizing Time Optimal Control Systems," Journal 
of Math. Analysis and Applications , vol. 1, pp. i|8ii-ii93, 1963* For the research 
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study reported in this paper, the proposed method was adapted for analog 
computer solution, and the regulator problem, x(t^) * 0, for three 2-dimension, 
time invariant systems, a neutrally stable system, a stable system, and a 
lightly damped system, as well as one 3-dimension, time invariant system were 
investigated with some analytic work developed to aid in the problem analysis. 
The computational procedure for the above mentioned systems as adapted 
for analog computer solution can be briefly outlined as follows: given a 

control system as described by equations (i-l) and (1-2), define 
(1-3) J = X T (t Q ,t) yj 
then 



d-t) 



-A T (t) J, j(0) - TJ 



Take some initial guess for yj such that Y) • x(0) 4 0, solve this set 
of equations to obtain j“(t) from which, using equation (l-2) 

(1-5) u- sgn (BT(t) J(t)) 

Now form the function 

(1-6) f(t ^yj) = Yi . (t( t l9 yj) + 3c(0)) 

where by definition 



(1-7) 



■*«, , -i, V 



X(t Q ,s) B(s) sgn (B T (s) X T (t 0 ,s)>^) ds 



so that ^ 

(1-8) f(t x ,^) = I ||^ T X(t Q ,s) B(s) | ds + =7-^(0) 

Determine the stopping time, t^, from the zero crossing of f(t^, Yj) for 

t = t-^ when f(t-^,>^) = 0, Now, 0 ^ t-^ ^ t-j 0 , where the time t-j° is 

the minimum time to drive the system to its null state. Obtain the next 

trial value of rj from the correction vector 

(1-9) 4vj - (*(t ,^f) +x(0))At 

where z“(t^, Yj ) is obtained by solving 

(1-10) y = A(t) y + B(t) u , y(0) =0 
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up to the time t^, and then solving backward in time 
( 1 - 11 ) w - A(t ) w , "w( t^) * y(tj) 

which yields 

(1-12) w(0) « 'z(t 1 ) 

As long as there is a correction vector, a new trial vector Y) may be 
used to make t^ approach t^° and YJ to approach Y) ° where >7 0 
is the constant vector which yields the optimum control vector u°. 

The outlined procedure can be adapted to the more general problem, 
however, the computer capacity required would necessarily be considerably 
extended. For this investigation of the general suitability of the 
computational method involved, it was decided to first investigate only 
the more simple types of problems already mentioned. 



3 



II. STATEMENT OF THE PROBLEM 



Given a control system, assumed to be normal, with constant 
coefficients and a single input, l) determine the optimum control law 
on the differential analyzer using an iterative technique proposed 
by Lucien Neustadt, 2) investigate the influence of the initial 
choice of the constant vector rj , 3) investigate the effect of the 
size of the sampling interval AT , and 1*) determine whether or not, 
for finite values of AT , the constant vector Yj might leave the 
domain of Y7 , i. e., whether or not • x(0) > 0 for any iterative 
value of Yj . 
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III. DISCUSSION 



The solution to the classical regulator problem as researched for 
this paper follows with minor exceptions the method presented by Lucien 
Neustadt ( Journal of Mathematical Analysis and Applications , vol. 1, 
pp JU81i— ii93 ) • Given a control system such that 
(III-l) x * A(t) x(t) + B(t ) u(t) 

where: A(t) represents an ri x n matrix, continuous function of time 

B(t) represents an n x r matrix, continuous function of time 
x(t) represents the n-dimension state vector 
u(t) represents the r-dimension control vector 
what is the optimum control, u°(t), which will drive the system from 
any arbitrary state, x(tg), at time tQ to some desired final state, 
x(t-^), at time t^ in the minimum time. It is assumed that the admiss- 
ible control is magnitude limited, |uj ^ 1, and that the system is a 
normal system, A system is normal if any component of B^(t) X^(tQ,t)Y] = 0 
on any finite interval of time, O^t^oQ , implies ^ = 0, where Tf 
is an n-dimension constant vector and X(t,tQ) is the n x n transition 
matrix for 

(III-2) 1c = A(t ) x(t) 

satisfying the equations 

(III-3) X(t,t 0 ) = A(t) X(t,t 0 ) 

(ni-lt) x(t 0 ,t Q ) = 1 

For such a normal system the optimal control is unique and is bang -bang 
with u^(t) * ±1, where the sign is given by 
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(III-5) u(t) = sgn (B T (t) X T (t 0 ,t)^) 

The solution for the system of equations described in ( III-l) is given by 



where C(t) has the properties: 

1 ) compact and convex 

2) C(t) D C(t^) whenever t > t^ 

3) C(t) grows continuously with time 

For t = t-^ 0 , -x(0) 6 . C(t-^°) and is a boundary point of C(t^°). Since 
C(t^°) is a convex set there is at least one row vector, such that 



so tnat ior x^t ; = u 




Define the subset GQ t ^ : 

(m-8) 



(III-7) 





(III-9) - rj . x( 0 ) > "9 *60 , for all co £ C^ 0 ) 

Thus "rj • cS takes on its maximum value when cO = -oc(0), see Figure III-l 
for a two-dimension example. 




\ 



Compact Set C(t °) 
Figure III-l 
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Now 00 e C^ 0 ), SO 



that 

(V 



(III-10) 

and 

(III-ll) 



X(t Q ,s) B(s) li(s) ds , 



uj ^ 1 



OJ = 



yj T X(t Q ,s) B(s) u(s) ds , | u- i | ^ 1 



so to maximize • to , so that ^ • cO = - vj . x(0), let 
(III-12) ^ «= sgn (BT(t) xT(t 0 ,t)yf) i 

Now define 



(III-13) 



z(t 1 ,'7 ) 



X(t Q ,s) B(s) sgn (B t (s) xT(t Q ,s)^) ds 



so that "z(t^,aj) is on the boundary of C(t^), see Figure III-2* According 
to equation (III-9) 



(IIX-lU) ^ • ’a(t 1 ,vf) > ■*? • t > for all J 6 C^) 5 

J t z'(t ]L , yj ) and the function y] ♦ ^(t^ vj* ) is a non-negative, non- 
decr easing monotonic function of time. 




Compact Set C(t^) 

Figure III-2 

Assuming the origin can be reached in some time, t °, so that -*x(0) 
lies on the boundary of the set C(t^°), there is a convex set of vectors, 
Hq, such that if yf £ Hq, - • x(0) maximizes the function q • tO for 

£ c ( t i°) and -x(0) * ^(t-j, 0 , ?))• Consider the function 
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(III-15) f(-t, ">7 ;x(0)) « yj • (z(t ,r[) + ~x(0)) 

Now by definition, ■ 0, so restrict Yj such that 

(III-16) f(t Q ,^ ;x(0)) < 0 

If £ H^, fCt^ 0 , 1 ^ j3c( 0 ) ) > 0, hence for some time t^, 0 4 t^°, 

(HI-17) fCt^vj ;x(0)) = 0 

Let 

(III-18) FCr? ,x(0)) = t x 

so that 

(III-19) f(F(vJ“ ,x(0)), v? jx(0)) * 0 

Since f (F(>? ,"x(0) ), yj ;x(0)) is continuous in its arguments, F(yj,x(0)) 
must be continuous in “rj so that if ^ Hq, F(*Vj" ,lc(0)) ^ t^° and if 

F(vy ,x(0)) = t^°. Thus any vector yj such that V] • *x(0) 0, 

which maximizes the time t^, for which yj • (z(t-^ 9 rj) + x(0)) - 0 
where is given by equation (III-13) may be used in the optimum 

control function. Conversely, if v] defines the optimal control function 
it maximizes the time. 

To solve for use the method of steepest ascent. Assume 

(III-20) rj = rj(V) 

then 

(HI-21) d^/dlT = k VF(y],x(0)) 

where k > 0 is some constant or function of IT and vj . In order to 
obtain a valid solution: 

1) F(V| ,x(0)) must have a maximum for "yj° 

2) the solution v^(tr) never leaves the domain of F(T],x(0)), and 

3) the partial derivatives exist, where are the 

components of yj. 

Proceeding formally, see Neustadt's work for statements of proof, 

from equation (III-19), given a specified initial state, 3c(t ), 

0 
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where F(>?,x(0)) = is given implicitly by (III— 19 ) • Then 
(III-23) V? 






but 

(III-21*) 

and 

(111-25) 

so that 
(II 1-26) 



C>f 



£= = (zftpY) ) + x(0) ) 



ii £ f 



Ib T (s) X T (t 0 ,8)-i] ds) 



VF = yj) + x(0))/ 

From equations (III-18) and (III-20), 

(III-27) F«F(^(f)) 



B T (t) X T (t Q ,t) 7f 



hence 

(111-28) 

let 

(I 11-29) 

then 



if _ ^ _ V7 rr 



«3l = 



(1 11-30) -(z‘(t 1 , Y] ) + x(0) ) 

so that finally 

(III-31) ^ -(zft^ + x(0)) df, 1) 2 - v^ l = dv^ 

To compute this corrected vector, , it is necessary to precompute 

—V. . 

z(t^, Y]^) for which in turn it is necessary to precompute the tune t^. 

Thus it is seen that "z(t-^, is not an instantaneous function of time 
and a finite sampling interval, AT' , is required, hence equation (III-31) 
becomes 

(III-32) \ ^ - (z(t 1 ,N i j 1 ) + x(0)) AT 

The time t^ was defined by equation (III-17) and "z(t^, ) is given by 

equation (I 11-13) so that we have 



B T (t) X T (t Q ,t)>)|| VI 
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(III-33) 



= J X(t Q ,s) B(s) sgn (B t (s) X T (t o ,s)>^) ds 



or using the steering command given by (III-12) 

(III-3U) ^ . i(t 1} VJ 1 ) . f h | B T (s) X T (t 0 ,s)v^ i || ds 



Now let 

( 111 - 35 ) 

then 

(111-36) 



J(t) = x T (t 0 ,t )^ 1 



I - -A T (t) J(t) , J(t 0 ) - ^ 

Solving equation ( 111-36) then gives *j(t) from which, together with 
equations (III-15), (III-17) and (III-31;), it becomes possible to compute 
the time t^. See Figure III-3 for a block diagram description of the 
computing technique. As revealed in Figure III-3, the control signal 
given by equation (III-12) is easily developed at the same time. 

Now with the time t^ known it becomes possible to compute "z(t^, Yj )• 
Consider 

(III-37) y = A(t ) y(t ) + B(t) u(t) , y(0) = 0 
which has the solution 



(III-38) 



y(t) = x(t,t Q ) 



X(t ,s ) B(s) "u(s) ds 



which combined with equations (III-12) and (III-13) yields 
(III-39) y(t x ) = X(t 1# t 0 ) t( t^) 
or alternately 

(111-1*0) t(t 1 ,v] 1 ) = X" 1 (t 1 ,t Q ) "yC^) 

where X”^(t^,t^) is the transition matrix for 
(111-1*1) "w = A(t) w(t) , 'w(t ]L ) = y(t 1 ) 
run backwards in time so that 

(111-242) V(t 0 ) - X _1 (t 1 ,t 0 ) yf^) 

= zttpV^) 

See Figure III-l* for a computational block diagram for the" function "z(t^,v^) . 
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IV. COMPUTER DESCRIPTION AND CIRCUITRY 



The analog computer used in the research of this paper was the 
Michigan twenty amplifier computer. The amplifiers had a gain of 
100,000 to 1. Amplifiers were not drift stabilized, however, accurate 
balance was accomplished through use of a 0.1 volt full scale deflection 
voltmeter. Resistors were matched to provide an accuracy of 0.1$ and 
capacitors were variable and adjusted to provide matched outputs. 
Twenty-three ten turn helipots plus a three gang helipot were set 
against a null heliDot to provide accurate potentiometer settings. 
Individual problems were set up on removable patch boards. 

In addition to the basic computer, a four place digital volt- 
meter was used to read out voltages. Automatic hold switching was 
accomplished by using relays in external circuitry. 

The circuits presented in this section were developed for the 
general two-dimension, single input system: 

(IV-l) x = A(t) "x(t) + b(t) u(t) , | u | ^ 1 

The three-dimension, single input system circuitry requires little 
modification and will not be presented in this section. For multiple 
input systems, a normalization scheme other than the one used in this 
work would be required. The specific circuitry for each problem 
analyzed is presented under the appropriate sub-section in Section V. 

Given the system described by equation (IV-l), the adjoint set 
of equations is given by 
(IV-2) 3 = -A T (t)j, 3(0) 
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where 



a ll 

a 12 



a 21 

a 22 



( IV-3 ) A T (t) = 

The circuits required for the solution of this set of equations are 



simple integration circuits as presented in Figure IV-1, where 

l*uJ 



P, = R n C 



a u F n 



~<t 



and 



r 2 c F 1 ■ 


1 a 2ll 




o< T 

0 1 


°^t 


R G 


| a 12| 


. 'J 1 


3 F 
* 2 


~^JT 


°<t 




l a 22l 




t 





Pj_ * potentiometer setting 

R^ = input resistance, ohms x 10^ 

C F = capacitance, farads x ICT^ 

cxi - physical quantity 

e e voltage representation 

o' - Veal time 
^ "^computer time 



B T (t) X T (t 0 ,t)lf =B T (t)j 



(IV-!*) 

or for the single input system 
(IV-5) 

a simple summing circuit as shown in Figure IV-2 is required for the 
solution of the set of equations given by equation (IV-£), where 



B T (t) X T (t o ,t) ^ 



o< 



7i 



B'I(t) X'I(t 0 ,t)Tf 



*6 

-HT 



2 I 



c*L 



J 2 



: B T (t) X T (t 0 ,tH 



and 



Rp. ■ feedback resistance, ohms x 10^ 



lli 





FIGURE 22-1 CIRCUIT'S FOR ADTOIVT SET OF EQUATIONS 




Figure BT-E 



CIReuir FOR SOLUTION OF B T (0 X T (t./0 s ^ T J 



For the single input system 

(IV-6) || B T (t ) X T (t 0 ,t) >l|| = | B T (t) X T (t 0 ,t) VJ | 

hence the absolute value circuit presented in Figure IV-3 was used 
to normalize the scalar quantity. To normalize a vector quantity a 
different scheme would be required. The absolute value circuit is an 
extremely accurate circuit. 

Now then 

(IV-7) Yj • t - 



B T (s) X T (t 0 ,t) | ds 



Solution of this equation requires a simple integrating amplifier, 

see Figure IV-h, where 

p . p P °<B T (t) xT(t n ,t)^ 

7 7 f 3 cx y) • i <* t 

By putting the initial condition equal to + • x(0) on the amplifier 

the function f (rj /z$x(0) ) can be obtained directly. To obtain the 
stopping time, run the system until 

f (vj ,^;x( 0)) * 0 

at which time switch to hold. Automatic hold switching was obtained 
through use of a relay as depicted in Figure IV-£. To insure accurate 
switching, the voltages were read out to four decimal places, an 
adjustable input was supplied to the amplifier as shown in Figure IV-5 
and stopping times which satisfied f ( 7} /z$5c(0)) = ±0.000£ volts were 
accepted. 

Time was generated on a simple integrating amplifier with constant 
input. When the computer was run in reverse time the input voltage 
was reversed in polarity and the automatic hold circuit was switched 
to the output of the time generating amplifier. Again the system 
was automatically switched to hold at time t * 0. 
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FIGURE 17-3 ABSOLUTE VALUE tlttlMT 
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The control signal was produced by taking the output of amplifier 
number 3, Figure IV-2, and feeding that into the bang -bang circuit 
presented in Figure IV-6. The potentiometer is adjusted to produce 
the voltage representation of the magnitude of the control signal. 

The voltage potential of two common outputs will be that of the most 
positive of the individual amplifiers. At point A the voltage will be 
the most positive value of 1) + saturation, 2) - saturation or 3) u = -1 
which automatically eliminates 2) as a possibility, since if 
B T (t) X T (t Q ,t)>7 > 0, -1 > - saturation and if B T (t) X T (t 0 ,t) vf < 0, 

+ saturation > -1 > - saturation. Similarly at point B the voltage 

represents u - +1 for B^(t) X^(tQ,t) >7 > 0, or u = -1 for B T (t) X^(tQ,t) V7 < 0. 

This is an extremely accurate circuit for which the voltage representation 

for u was set accurately to four significant figures. The switching 

time was determined to be less than 0.1 milliseconds for an input 

changing at the rate of £0 volts/sec which was chosen as a representative 

rate of change of the voltage representation of B^(t) X^(tQ,t) v? • 

The circuits used to compute z(t , yj ) were simple integrator 
circuits. First y(t^) was computed, see Figure IV-7, from 
(IV-8) y « A(t) y + b(t) u , y(0) = 0 

and then z^t^,^) was computed by integrating 
(IV-9) w - A(t ) w , w^) * yit^) 
backward in time, see Figure IV-8, to yield 
( IV- 10 ) = w(0) 

To perform the backward integration, it was only necessary to reverse 
the polarity of all the inputs to the integrators used in the 
generation of y(t^) and remove the control signal input, all of 
which was accomplished by manual switching while in the hold condition. 
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FIGURE JS.-G BANG - BANG CIRCUIT 
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FIGURE ET-7 CIRCUITS FOR TH£ GENERATION OF* /*$ (t) 
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FIGURE DZ - 8 CIRCUITS FOR THE GENERATION OF Jir(t) 
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In Figure IV- 7 



P 8 = R 8 C F. °<t 

k 



P„ = R„C 



a 12| 



Z2_ 



9 9 F li ^ 7 ! 

I ^1 1 u 

P 10 “ R 10 C F j< =< y x ~t 



P_, * R,tC_ ^ 

11 11 Fr ^ Qx £ 

p J 2 



a 2l| <?<?-, 

y , 



■221 



P 12 ° R 12 C F 5 c< t 



| bp [ <X „ 



P 13 ’ R 13V <=<y 2 (X t 



and in Figure IV-8, 



a lll 



P lJi “ R Ui C F 6 o<T t 



P„„ « R,„C. 



a n | °< w „ 



« 15 F 6 



a 2 l| C><w, 

P., = R 1(C C_ I I 1 

16 16 F ? ^ cx t 



P 17 * Rn 7 Ctj 



22 1 
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V. 



SPECIFIC PROBLEM SOLUTIONS 



Selected two-dimension as well as‘ one three-dimension system were 
investigated for the research of this paper, ^he two-dimension systems 
investigated were: 

1) a neutrally stable system, x = ® u 

2) 



u ^ 1 



Stable system, ~x = ll ^ + ll^j U 9 ^ u 

L -v~ ^ f n il f n" 



3) a lightly damped system 
an 

h) 



. r ° 1 

’ [-1 -.2 



6 1 
u 



and the three-dimension system investigated was: 

0 1 Ol (o' 

0 0 1 x + \ 0 } u 

0 0 0 



x = 



u 



Plant 1) was chosen to provide a System on which a rather thorough 
analytic solution could be obtained in an effort to gain some insight 
into the problem solution. Plant 2) was chosen to investigate the 
effect of running a stable plant in reverse time, that is a stable 
and an unstable combination. Plant 3) was chosen as a representative 
two-dimension system for which solution might be desired. The three- 
dimension plant U) was selected to provide a simple extension of the 
computational procedure to higher order systems. 
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Plant 1) 

The analytic solution to the control problem synthesis parallels 



the analog computer solution. Given 



(v-1) 

define 
( V-2 ) 
where 
(V-3) 



0 1 
0 0 



0 0 

-1 0 



?(°) = 



$ 1 



^11 
'n 12 



3 = 



23 



X(t 0 ,t 0 ) = I 



then 



(V-U) B T (t) X T (t 0 ,t) yj ± - b T J - -i7 u t + ^ 12 

from which 

(V-5) u = sgn (-vj u t + yj ) 

and |'t i 

(V-6) • t - ) | -^ n s + Trj | ds 

t O 

Referring to Figure V-l, for and the same sign 

(v-6 a) h 12 | t-|T7 n )^?, t < v Jl2 /r }ll 

2 

(V-6 B) % -t - -I YJ 12 | t + |TT} U 1 , t 4 Y| l2 /Tf ln 
and for and Y} 12 of opposite signs 

2 

(v-6 c) T q 1 - t - 1^1 t + |y) u | , t > 0 

Now 

(V-7) f(^ 1 ,t;? 0 ) - 0 

where 'xq = 'x(O) - ^01^ , defines the time t^ implicitly. From (V-6 A) 

(V-7 A) t = ^1^ ~ ^ ^12' 2 * 2 '^lJ ^ 11 X 01 +> ?12 x 02^ 2 

1 IV 

from (V-6 B) 

(V-7 B) t - lV| * ~ 2 ^lll (> hl X 01 + ^12*02 + I ^lgj-^lul ^ 

1 111 

and from (V-6 C) i 

(v-7 C) t = ~ I ^ 12 1 + ^12 1 ~ 2 ^lJ ^n X 01 * ^12*02 *'^12^711^- 

1 I Rll I 

By definition 

(V-8) '^ t 1 ,T7 1 ) = j X(t Q ,s) B(s) sgn (B T (s) X T (t Q ,s) Ijfj) ds 

but 0 

(V-9) X(0,s) - ^ “®j 

S ° /' t 

(V-10) " Jo f-'l sgn (-Y^ 3 + Y) 12 ) ds 

so that for t ^ ^12^ 11* and ^11 2111(5 ^)]_2 of the sanie sisn 
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(V-10A) z(t i ,^ i ) = 



sgn (Y) 12 ) 
t x sgn (>] n ) 

for t-^ ^ ^12^ Y ^11 , and ^11 and ^12 sairie sd & n 

(¥-10 B) -z(t 1 ,^ 1 ) = /sjyo - kf'l.i/n,,) S<jVi (%) t '\ Sj* <\Y 



and for Y] and Y) of opposite signs 

'11 n ^ 



j r '"' 59Yi(p M ) t* 

c f W*,, 8 9 «(y, m ) _t sjxrriJ J 



(V-10 c) 

Now 
(V-ll) 
and 
(V-12) 
so that 

(V-13) 



'12 

d 



(t*)/2 sgn (1? n ) 
sgn (y) i;l ) 



W = -(«(t 1 ,^ 1 ) + * 0 } 
A \- 



\ + i m y ± - ^vV + x o ) ** 

With this computational procedure, successive values of Y} were obtained 
for AT = 1, x Q = [ 1} . The data are presented in Table V-l and a plot 
°f Yj^ versus is presented in Figure V-2. From equations ( V— 7 ) it 

is apparent that t^ is a function of only one variable. To gain a 
better understanding of the function F(^,Xq) on which the method of 
steepest ascent was used, it was decided to find F(P ) where P - Y^/y)^. 
Then for ^ 4 V] 12 A}u, ' )r ) 11 and ^22 the same si S n 
(V-U* A) tj_ - P - ( f> 2 - 2(x Q1 + P x Q2 ))i 

for t l ^ Y) 12 /^ 11 ’ ^11 and ^12 the Same Sign 
(V-ll, B) tj_ = P + (-/= 2 + 2(x q1 + P x 02 ))i 

and for Y) ^ and of opposite signs 

(V-ll* C) tj_ - P ♦ (p 2 ♦ 2(x 01 + p x 02 ))i 
F ^l»x 0 ) was plotted versus in Figure V-3. 

The effect of the size of A't' was investigated by linearizing equation 
(V-13) about T^° = ^ i which was obtained from Figure V*-3. Since Yj^ 
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TAQL? I l-X ArtALYTlC DATA FOR 





JUU 



flu>H 


>?. 




k 






M '/it 


d Ht\ 


1 **> 


| 4n. 


At- 


i.o 


















1 




. S' 


0 


0 


0 


-f.O 


O 


- 1.0 


0 


Z 




.S' 


1.0 


-S' 


f/.O 


-s' 


-/.o 


-S' 


-/.o 


3. 


-/.s 


-•S' 


/.706 


-/.jv/ 


+/.OWD 


■*.MZ 


-/.O+O 


i.3Hl 


-/.O+O 


V 


-US'S 


-LS'iO 


1.3/0 


i./zz 


-.360 


-t./jz 


+. 8 S 0 


-1./3Z 


+.3 60 


5 


-2.190 


-.690 


!.S9t 


-/. 3 VD 


+/.089 


+.3+0 


-1.039 


+.3+0 


-1.037 


6 


-l.lSo 


-/.779 


1.99 Z 


-US7 


+.168 


4.JS 7 


-/it 


+./S7 


-.163 


7 


-1.192 


-1.917 


1.99/ 


-.306 


-o7S 


-J9Y 


+.D7S 


-19+ 


+.07S 


8 


~/. 1 B 7 


-1.371 


1.997 


-l./o& 


+./// 


+.IOS 


-/// 


+./OS 


-. 1/1 




- 1.8 81 


-t.993 


1.99 S 


-.888 


-.in 


-. HZ 


+./// 


-.nz 


+./// 


/0 


-/.ffV 


- 1 . 8 ft 


/.Q96 


-/./ol 


+./10 


+./09 


-JLO 


+./09 


-JlO 


// 


-/.0W 


~/.31Z 


1.117 


-.376 


-Ml 


-lz + 


+-J/7 


~.IZ + 


+,i/7 


/2 


-2.609 


-1.3TS 


2.000 


'llll 


+.13 1 


+./tl 


-J3Z+ 


f./t? 








at: c 


o.z is 


















i 


-J.0 


0 


i.+H 


- 1.0 


tint 


O 


-/.V/Y 


0 


-631 


l 


-1.6 


-.43/ 


/.83? 


-1.W06 


9-.17 7 


+.+06 


-.717 


+./SZ 


-.11* 


3 


-.3+3 


-.«zj 


/.?»? 


-/. 06 J 


+.061 


+.067 


-.069 


+.OZI 


-ozz 


H 


-.317 


-8¥/ 


Z. 00 O 


-.967 


-.O+Z 


-•0+3 


i.OHZ 


0/6 


+.0/6 


s 


-SH3 


-.82? 


z .000 


-t.oSi 


+. 0 / Y 


+.03 + 


-o/Y 


+.0/3 


- 00 S 


6 


-.330 


-.«jy 


2.000 


-.190 


-.608 


- O/o 


+. 00 B 


-.oof 


+.OOZ 
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TABLE X-X (tONTtMLLZD) 



RUH 


* 




t, 


4,&> 


oM***i\ 




*1, 


Ay %. 


AX ” o.zi 


1 


- t.o 


o 


/.v/v 


-t.o 


U V/V 


O 


-zv/y 


o 


-3SJ 


z 


-t.O 


-.333 


A 722 


-I.3S7 


+1.0/4 


+.3S7 


-1.016 


+.089 


-.249 


3 


-.911 


-.407 


/.«/ 


-t.246 


4.S67 


4.394 


-.547 


4.0 99 


-.192. 


9 


-.8/2. 


"7 Vf 


mv 


-Z./J8 


4.1*3 


+J38 


-.198 


4.039 


-.03 7 


5 


.i 

CO 


-.786 


1.999 


-473 


-0Z3 


-0lS 


+.OZJ 


-.004 


+.004 




At s O.t 


1 


OAfA 


F0f 


At ’ 


to 












z 


0AT4 


FOR 


At?* 


f.0 












a 


-/.Soo 


-Soo 


/.7a 7 


- 1.396 


uni 


4. Hi 


- 1.191 


+.03S 


-.119 


9 


-LHiS 


-AH 


/.7« 




4.931 


4.284 


-.931 


4.039 


- 093 


5 


‘HU 


-.707 


t.tio 


- t.tol 


*.izt 


4. *09 


-U8 


4.091 


-.063 


6 


-/. 3tS 


-7 90 


/.»« 


-tHU 


4.7 28 


4.944 


-7l8 


4.093 


-.073 


7 


-t.391 


-.843 


t.403 


-/Hot 


4.6/9 


4.901 


-Mt 


4.090 


-.OU 


8 




-MS 


Ml 


-1.284 


4.3U 


4.384 


~.St3 


4.039 


-.031 




-t.ui 


-417 


I.ISt 


-t.zio 


4.90% 


4.210 


-.90% 


4.031 


-.09/ 


to 


-/.IZI 


-1.018 


HU 


-Kill 


4.313 


4.1V 


-.3lS 


+.027 


-.031 


it 


-t.ioy 


-t.oSo 


/.18S 


-l.tol 


4.1*1 


4.Z09 


-.191 


4.0Z1 


-.oil 


n 


-1.183 


-1.019 


t.99l 


-1.163 


4.180 


4.143 


-.180 


4.016 


-.01 & 



28 




FI6U.RC 3T-a PLOT OF Vl it y/CKSUS *Ju 




29 







- O 

II 

I * 0 
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FIGULR0 tf-J FCrf,^) VStSiiS e FC>K 



and < 0 and usin § equations (V-10 B) and (V-13) 

(v-13 b) >j i+1 = / - L - V^/ 1 -fW^,) 1 *■ a (Vo, * tiO 

define " + “ + 7 : 0U/*|J l +• afo, + >Mw + < t >ZT + \J 

(V-15) ' 



<r 



"Vi+1 " \ " <P 



d *h.+1 

1 



H -V ‘‘I 



2 /if 

(i - 2 Ar ) 



.] 



then for x Q =^q] »^° " {llj 

(V-16) f » A dij 

where 

(V-17) A - (1 - 2 2 At r) 

The eigenvalues of A are 

(v-i8) >1 - i, i - hAtr 

For stability, lA.I < 1, hence for the given plant with lc Q = J 
Atr < 0.50 gives stability. As noted in Table V-l and Figure V-2 
with ATT * 1 there is a slight divergence in the successive values of 

— ^ — V Q 

Y) as Y| was approached. 

Using the circuitry presented in Section IV as combined in Figure V-b, 
the computational procedure was investigated for values of At' * 1.0, 

0.375, 0.25 and 0.1 with Xq = QJ- . The data for the iterations are 
presented in Table V-II and are plotted in Figures V-2, V-5, V-6 and V-7 
respectively with the corresponding analytically computed iterations. 

As expected ATT = 0.25, A = 1 and 0, yielded the most positive convergence 
to yj°, AX - 0.375 gave rapid convergence with some oscillations, 

ATT « 0.1 gave slow convergence and At = 1.0 gave large oscillations. 

The hill climbing pattern associated with each value of ATT is presented 
in Figure V-8. It was noted in Table V-II that as "x — ■> 0, the correction 
vector Cz(t^,> 7 ^) +'Xq)AT' — > 0 independent of the choice of AT: • The 
plots of X£(t) versus x^(t) are presented in Figures V-9, V-10, V-ll and 
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FlGULRt TT-W SCKfAIATlC DIAGRAM FOR 01 \/t + \°Cao 



RUN 


*1 




t, 




/*A)I 


3.M 






[*Sfe 






P 














X 


© 


.5 


o 


3 


o 


O 


O 


-i.o 


0 


-/.0 


27 


oO 


z 


-/.o 


5 


1.7 Oi 




9/.007 


”.5oo 


Hoo7 


- 5oo 


-Aoo7 


-.600 


“l.ool 


-.Soo 


$ 


“Moo 


~*>7 


l.7o$ 


./.vyf 


4t.es ( 


*/3S/ 


t/.OS 9 


V.35/ 


-/.039 


4.3SI 


-/.oJ9 


4.368 


H 


“/./Jt 


-/. 6 Y 6 


118(0 


-.V/9 


-954 


4.17$ 


-9Y2 


-/.476 


9.9Y2 


“i.ZIS 


MV 2 


H.ity 


s 


-2.Y/V 


-. 6 oY 


UYf 


f/jra/ 


9/./JJ 


•us* 


3/.Z37 


4.Z87 


-/./J 7 


4.2 89 


-/♦/37 


4. ISO 


6 


-Z./25 


-/. 7 V/ 


Hit 


9.346 


*.322 


-LU$ 


*.326 


4. Hi 


-.326 


*♦ 2/6 


-.326 


9.8/9 


7 


-W 1 


-2.043 


2 .ooo 


~.ti$ 


-.180 


i 

2 


-./78 


-.14$ 


*./78 


-./?5 


f./78 


4l.o8o 


3 


- 2 ./oy 


-/.tg5 


Ml 


+./9Y 


bISS 


-/./66 


4jS8 


4. US 


-./88 


*./66 


-./ad 


t.tu 


9 


-/.931 




/.88J 




“J6V 


-.82* 


-/5J 


-m 


*./5J 


“I7Z 


4.16$ 


4/.0H 


/O 


'2.//I 


-/.no 


2.oo3 


4.ni 


t.uz 


-/./J9 


f.662 


4.137 


-./42 


4.(37 


— . /6 2 


4.910 


// 


-/.972 


-2.082 


2.009 


-.11 4 


-.us 


-.84/ 


“.US 


-./J9 


4.US 


-.13*7 


4. IZS 


H.oSS 


/2 


-2.//I 


-/•9$7 


ton 


9./J/ 


hno 


-/.i/5 


4./Z 9 


4.//S 


-/29 


4. ns 


“IZl 


4.9 26 


Ot'O 3W 


1 


-/.o 


O 


/. YOf 


t/.ffJ 


t/MOt 


-.959 


HMlo 


“.01/ 


-I.4Z0 


-.004 


-.SJ3 


O 


z 


-1-ooi 


~.$J$ 


/.«v/ 


41.01/ 


4.170 


-/.39/ 


4.770 


4.341 


-.770 


4.147 


-.288 


9.6J/ 


J 


~.ts7 


-n\ 


/.99J 


+.07S 


4.07 / 


-/o65 


4-OV 


t.oSS 


-.07/ 


4.oll 


-.027 


A 916 


H 


-•35 


-iva 


/.99J 


-039 


-.0*/ 


-.95o 


-.OYJ 


-.050 


4.046 


“.0(7 


4.0 li 


F/.OI? 


s 


-.WV 


-.9J2 


/.999 


4 . 0 HS 


*OVJ 


-/.otf 


9.0VJ 


4.0ZS 




+.611 


“.0/S 




6 


-.ftv/ 


-*v« 


2.000 


-.043 


-043 


-.96? 


~.blS 


“031 


4.01$ 


-.0/6 


4.009 
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TA©l£“ r-IT (cotiTlNULCo) 



RUV 




•*» 


rr 


•*A) 1 J I <3. (O 




^rl 


U,| 






At'* 


C 2 A 
























A 


-/.0 


O 


/.vo? 


4/.9?3 


4/.V08 


-.98? 


+1.420 


-on 


- 1.410 


-.002 


-.355 


O 


e 


7.002 


-.355 


/.72/ 


4/.JZ8 


+ 1 . 00 S 


-;.359 


tl.ooS 


t.359 


- 1.00 5 


4.087 


-.252 


4.35V 




-.9/5 


-.407 


/.9/V 


+.730 


+.976 


7.385 


♦.578 


4.385 


-.578 


4.09/ 


-./V5 


4.47V 


v 


-.82V 




/.w 


+.161 


+.13 8 


-tw 


4./5B 


+.147 


-./58 


4.0J7 


-.03 9 


+.912. 


5 


-.7*7 


-.79/ 


2.ooo 


-.023 


-OZJ 


-.949 


-.025 


-.031 


4.025 


-OIL 


4.009 


+ 1.007 


LZ * O.i 












1 


o/\r/ 


FOR 




> 


















2L 


DATA 


FOR 


Af e i./ 


> 


















3 


-l.Soo 


-.5o 7 


1 -ToS 


//.Wf 




7.35/ 


4/.OJ9 


*3S1 


-7.039 


4.035 


-.104 


4.338 


V 


-1.463 


-.4// 


/.74S 


H.Z58 


+.933 


7.375 


4.933 


+.378 


-.933 


4.0JS 


-.09J 


+.41 7 


5 


-1.417 


-70V 


/.8Z9 


+/.//5 


4.8V& 


-1.407 


4.832 


4. VO 7 


-.832 


+.041 


-.083 


+.H9J 


6 


7.384 


-.787 


2853 


4.935 


4.725 


-7.406 


4.725 


4.Y05 


-.725 


+04/ 


-.073 


4.548 


7 


-I.34S 


-.840 


/. 9/2 


+.778 


4.4 zo 


-1.406 


+.6lo 


+.406 


-520 


+.04/ 


-.04 4 


+.640 


8 


-/.3*V 


-9/2. 


1.93/ 


4.4VJ 


+.522 


-1.178 


+.S21 


4.3 73 


~.5Z2 


+.038 


-052 


4.499 


9 


7.244 


-.94V 


/. 9VZ 


+.5o£ 


■/.Vi 7 


-1.130 


f.4l7 


4.3 JO 


-.V*7 


4.03J 


-.043 


+*.750 


/O 


7.233 


7.007 


/.?6? 


4.382 


4.337 


7282 


+.336 


4.282 


-JJ4 


+.oz8 


-.03V 


4.8/? 


// 


7.205 


-1.04 0 


/.982 


+.280 


4.255 


7.ZV7 


+.W7 


+ .147 


-.257 


+.OZS 


-.016 


4.84 2 


/e 


7./80 


<0 46 


M89 


+/B4 


+.179 


7/58 


4./80 


+.168 


-./80 


+.617 


- 0/8 


4.902 
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00 
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vS 

*0 
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FIGURE 
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FIGURE: 7 -11 PLOTS OF V£*S<AS /*, 



4i 




FIGURE TT-IE 
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FOR 



V-12 respectively to indicate the relative error in x(t^) for successive 
iterations of rj • The initial choice of J on the boundary of 

Yj^ . x( 0 ) ^ 0 was chosen for At” = 1 . 0 , for AT = 0.25 and 0.37£> 
vf i was chosen and for At = 0 . 1 , > 7 ^ was arbitrarily chosen. 

There was no apparent effect due to this limited choice of initial values 
of v? r 

The plant was further investigated for another given initial state, 
*°> ■ til . The plot of F( >*q) versus Z 3 is presented in Figure 
V-13. Again the difference equation (V-13) was linearized about 

( V t - (-0.668 At +1) 1.198 AT 

6 ' 1.11*2 AT (-2.01*7 Alt + 1) * 

so that A has the eigenvalues 
(V- 20 ) ^ = 1, 1 - 2.715 Atr 

Analytic and computer computations for AT = 0.1 and 1.0 were made and 
the data are presented in Table V-III and displayed in Figures V-ll* and 
V-15. The hill climbing patterns are presented in Figure V-l 6 . Again 
AT = 1.0 gave oscillations while AV = 0.1 gave slow convergence, and 
as before as x(t^) — ^ 0 , the correction vector, -Cztt^,^ ) + Xq) A'C 0 
independent of the choice of AT. The plot x^(t) versus x^(t) is 
presented in Figure V-17 for AT = 1.0 to display the relative error 
in "x(t^) for successive values of > 7 . 



Plant 2) 

The computer was set up as shown in Figure V-l8 using the basic 
circuits of Section IV for 



(V-21) 



-2 

0 



lul ^ 1 



Initial states of x(0) = , ^3.5j ’ and [ 3 o] were investigated 



1*3 






•ft 
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TABIC 2 -nr ANALYTIC AND COHPU,T£H DATA 



RUM 


7, 




t, 


^ ft.) 




3/0 




Sfc| 


L* 




1 ^ 


At * 


4.0 


Ammo 




















1 


-/.326 


-o< 2 Y 


/.wy 






-/.OVO 


9/. 330 


4.0Y0 


-.830 


4.040 


- 8 S 0 




2 


-A 266 


— .9/Y 


/.5"9J 






-767 


+./77 


-.233 


+.32J 


-.233 


+.323 




3 


-Z.5/9 


-.589 


/.597 






-/./ 2 / 


9.82/ 


+./Z 1 


-.32/ 


4./Z/ 


-.32/ 




H 


-/.J90 


-.9/0 


/. 6 // 






-.876 


+.Jo9 


-.119 


+./9/ 


— ./ 2 Y 


-./9/ 




3 


-2522 


-.7/9 


7.6/2 






-/.076 


4.668 


+.076 


- /68 


4.078 


~./ 6 d 




6 


-/.vw 


-.867 


/. 6/8 






-.933 


4.390 


-.067 


+.//0 


-.067 


4.//0 




7 


-/. s;/ 


-777 


/.6/9 








t.iol 


4.033 


-. /o? 


+.055 


-./*9 




8 


-MS 6 


-.866 


/.</? 






"9 VO 


+yo / 


-.060 


+.o?7 


-.060 


+.099 




9 


-iS/6 


-.787 


/.620 






-/.ov? 


4.380 


4.0Y2 


-.080 


4 . 0 V 2 


-.080 




JO 


-/.V7V 


-.867 


/.62o 






-.963 


4.WV 


-,03V 


+.AS 6 


-.OJY 


+.056 




n 


-/.So* 


-.fc// 


/.620 






-l.oiz. 


4.3Y2 


+. 022 - 


-.092 


4.022- 


-OY2. 




it. 


-/.VCfc 


“85*3 






















AlV * 1.0 


COnPiLTER 




















4 


-/.J26 


~.oui 


/.V60 


+/.zzy 


4.8/8 


-J.OSZ 


+ /.3J4 


-4.052 


-.836 


4.o3Z 


-.836 


4.0V7 


2 


-Z27J 


-.9/8 


/.6/9 


-.296 


-3jy 


-.732 


+.W5 


- 2 Y 8 


+.3J5" 


-. 2 Y 8 


4.333 


4.921 


3 




-.363 


I.AOl 


4.Y^ 


+.3V« 


-/./vy 


4.8Y9 


+. /vy 


-.339 


+./YY 


-35? 


4.370 


V 


-/.378 


-.942. 


/. <09 


-./96 


— ZlZ 


-. 8 vy 


4. 7 <6 


-./36 


+.ZJY 


-./5< 


4.23/ 


+.6 7/ 


5 


-(.fw 


-.668 


/.</f 


4.23/ 


4 . 2/0 


-/.086 


4.7/0 


4.066 


-. 2/0 


4.086 


-. 2/0 


4.YY9 


6 


-/.VV 8 


-898 


/.<33 


-.096 


— ,//3 


-.922* 


+.37J 


-.078 


+./Z 8 


-.678 


+.116 


4.62/ 


7 


-J.SZ3 


-.77/ 


/. 6 JJ 


+.//?9 


4./OY6 




<606 


4 . 0 S 2 


- /06 


4.052 


-JO 6 


+.306 


8 


~l Y7V 


-.876 


/.<32. 


-.0948 


-.0709 


-.96V 


v-.yzv 


-.036 


4.076 


-. 0 J 6 


+.076 


+ 59 V 


9 


-/.3/o 


-. 80 / 


/. 6 J 0 
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to determine the effect of computing solutions through different multiples 
of the system time constant. 

For x(0) = [°- 2 J , AT = 0.5 and 1.0 were used. The data are 
tabulated in Table V-IV, Plots of versus are presented in 

Figures V-19 and V-20, the x^(t) versus x-^(t) plots are presented in 
Figures V-21 and V-22, and the plots of yg(t) versus y-^(t) and W 2 (t) versus 
w^(t) are presented in Figures V-23 and V— 2I4. to indicate the range of 
y(t) and "w(t) in the computation of z(t^). As was the case with Plant 1), 

AtT = 1.0 gave some oscillations and relatively slow convergence, whereas 
AT =0.5 gave much more rapid convergence and eliminated much of the 
oscillations. Again it was noted that as "x(t^) — > 0, the correction 
vector -(z^t^,^) + "Xq)aT — > 0 independent of the choice of AT . 

For x(0) * , the choice of AT was dependent upon the 

criterion |Av^ + -jJ <[ |A^]^| • The data are tabulated in Table V-IV, 

the versus plot is presented in Figure V-26. Using this criterion for 

the selection AT rapid convergence with some oscillation was obtained. 

Again as the state vector ”x(t^) — > 0, the correction vector 
-(ztt^Y} ) +x q )AT > 0. 

For x“(0) = ^he choice of AT was dependent upon the criterion 

A^ + ]_ * 0* The data are presented in Table V-IV, the 

versus plot is presented in Figure V-27, the x^(t) versus x^(t) plot 

is presented in Figure V-28 and the plots of y 2 (t) versus y^(t) and w^(t) 
versus w^(t) are presented in Figure V-29. It is noted from Table V-IV 
that the additional requirement that V| • Xq 0 was used in the selection 
of AT, hence using the finite difference equation it does become 
possible for r | to leave the domain. In Figure V-29 it is noted that the 



52 



Rti* 






* 










1 








o? 4±2LO 










1 


-/.0 


O 


.7760 


+./V7? 


+.3637 


-./7/0 


+. 62 /? 


-.0070 


-62/7 




-.0090 


-.62/7 


l 


-/.007 


-.622 


.5750 


-.0275 


“2209 


-.0780 


-3760 


•JO CO 


+.3760 


io 


-JO co 


+ 3760 


3 


-4/// 


-.276 


.5?6Y 


t.OTZS 


r.2/7J 


-/ 6 oY 


V.3#t2 


-.os 76 


-.3882. 


l.o 


oJf 6 


-J W 


V 


-/./6 1 


- 4/V 


.6377 


-0/78 


-.077/ 


+.0678 


-1160 


-2o78 


+./760 


1.0 


-.2078 


+./760 


s 


-/.J» 


-.YJf 


.63// 


+.0372 


+.//7J 


-.237/ 


+.L/7V 


+.037/ 


-.1/7Y 


1.0 


+•.039/ 


-.2/77 


6 


-/.3/Y 


-.655 


. 6 VJ 0 


-ouo 


“0528 


-./v 77 


- 0?«2 


-032/ 


+.078 21 


1.0 


-.052/ 


+. 0 f« 2 > 


7 


-/.J 66 


"557 


.69*9 


f.O // 2 . 


+.0928 


-2o76 


+080S 


+.0094 


-.0805 


1.0 


+.0076 


-.0805 


8 


-/.JS7 


-.638 


.6723 


-.0076 


-.0276 


-/7VV 


- 0 V /8 


-.0236 


+.67/ 8 


1.0 


-.0256 


4.07/8 


? 


-/.38t 


-576 


.6Yo5 


t.OOSO 


V.O 203 


-,t/o 7 


+.0370 


+. 0 / 0 ? 


-03 fo 


1.0 


+. 0 / 0 ? 


-.0390 




*Llfo} at ; Q .5 


2 


-J .0 


O 


.Y?V5 


+./V76 


+.3863 


-/7o6 


+.6162 


-004*/ 


“6262 


0.5 


-.0077 


“3131 


e 


-l.ooS 


-.3/3 


.6267 


+.0378 


+./Z 8 S 


-.2676 


+2370 


+.0674 


-.23?o 


0.5 


+ OJ25 


-.U8S 




-MM 


-.VJ2 


.69VY 


-oo// 


-0060 


-./*?* 


-0076 


-.010$ 


f.OOfi 


a5 


-.oo5Y 


+0078 


V 


-M7$ 


-.727 


. 6 vVf 


tool/ 


+. 008 / 


-.Loot 


+.0/5J 


+.ooo7 


-o/55 


o.5" 


b.oooi 


-0057 


ar 


•Ml Y 


-.Wil 


,6Wf 


+.«*>? 


+. 002 ./ 


-:/fVo 


*.oo /6 


-.0060 


- 00/6 


as 


-.0030 


-oo# 8 












*■ iVi 


1 


~/.0 


0 


/. 28 * 




t.llss 


-3.6/2 


+2.706 


+.//Z. 


-2.706 


o.y 


+.056 


-/.35J 


2 


-.?vy 


-1.3S3 


/.5J3 


+.0777 


t.j(vV 


-9.i79 


+.772 


+207? 


-,7ft 


O.lS 


+i 26 ? 


-./?• 


3 


-.475 


-/.55/ 


/.i*l 


-.0694 


-.1170 


-/./7? 


-/.0/5 


-2.32/ 


+/.0/5 


0.25 


-.680 


+.256 


V 


-/.*« 


-/.Mi 


/.V87 


+.Z620 


/.3oo7 


-V.835 


+/.3J0 


9/.J3i 


-/. JJo 


0.25 


y.ijv 


-.j?* 


3 


-.?*/ 


-A627 


/..iv9 


+.0065 


9.M<« 


-J.670 


4.09*9 


+./70 


-ofof 


0 ./*i 


+. 0 Z/ 


-0// 


6 


-,?o 0 


-/.4J8 


/.5V7 


-00/3 


+.oo55 




V.ooY? 


t.ooY 


-.ooyf 














-A 

/IT.® 


ffl _ 




1 


-1.0 


0 


2. HI 


+.77Y 


+.870 


-J2.67 


+ZfV 


+ 2 . 6 ? 


-7.9Y 


,o«r 


Y./68 


-.Vf6 


e 


-wt 


-.Yt6 


Z.Z1L 


+.627 


f.77? 


-37.03 


+1ZS0 


+ 703 


-7.35 


.0/56 


+.//2 


-./2J 


j 


-.720 


-.6/7 


2.163 


+.i</ 


+.72/ 


-37.70 


+ 7.06 


+9-9/? 


-7.06 


.0/56 


y-./iV 


-//0 


i 


-546 


-727 


2.3oo 


■f.YfrV 


+. 6 Yf 


-YJ./f 


+6.60 


+/J./7 


-6.6o 


.0/56 


+.206 


-/03 


s 


-.340 


-.832. 


2.18/ 


+.3/6 


+.5o3 


-i/.l? 


+55? 


+ 2 /./? 


-i.-S9 


.0078 


+./6S 


-07J 


6 


-/Y5 


-vs 


z.V 8 o 


+.38/ 


+.067 


-Y8.77 


+2.33 


y/*9? 


-21.33 


.003? 


+.072 


-oo9 


7 


-J 2 J 


-.MY 


2.Y63 


-.027 


+.0/7 


-32.7? 


/.379 


+2.77 


*:i7Y 


,0020 


+.006 


-.00/ 


8 


-.153 


-•*/ 


21 V 6 3 


.ooo 


+.00 5 


'3/.V2 


+ .27? 


+ /.Y2. 


-27? 


.0020 


+.©o3 


-.00/ 



TAOL£ y-IT computed data FOR % - £-| | ®"J AJL 



55 



t 




54 



PLOT OF Y]^ v'tfflsus r \ POff 7t - [o J] ^ Cl) ^ ; e £ ‘^j ,AT*O.S 



>!.■* 




55 




56 








* r -lo 



57 




58 



I 




-/.* -u -4* “/.« -.1 -.7 -4 -* -.V -J -2 -I 



* 

.2 ^ 



.1 1 



FI6URF T- Z5 

PL or or ^£«»as 17, 

fo * 



OfTlliUM 



/ 



J-l 



■ jfc . 

'f * 



-i 

-.i 

-.2 

-.V 

-:S 

-.4 

-.7 

-f 

-1 

-/.« 

-a 

•u 

-/.j 

-*r 

-/.* 

-/.i 

-u 



60 




61 



- .» 



-.6 -5 


-.¥ 


-j 




r/ottRe 1 y* 


-27 




nor of ij At 


vm us 


Vm 


FOR ^’[' Z !]^ + [, 


- 




Lo-lJ v . 1 


'J 






-I 



opnmaN 



t 

'I, 




-.1 



-.2 



7* 






-.6 



y * 



-.7 



?S 



x « 
* 7* 



I 

I 

it 

/ 



?7 



62 



t 




63 



*- V ft 




64 



range of y(t) is much smaller than the range of w(t). Any noise in the 
system would have a more pronounced effect on y(t) resulting in errors 
in y(t^) and hence in'z'Ct^). This could possibly explain why 

) + x 0 )t^ 0 as x(t^) ^ 0* This problem, or source of 
probable error, could be minimized by rescaling the amplifiers used to 
compute y(t) and w(t) when switching from forward time to reverse time, 
or by using separate amplifiers to compute the two functions* 

The only analytic work performed for this plant was the computation 
of the time t and the results are presented in Table V-V. As can be 
seen from the table, there is good agreement in all cases between 
the analytic and computer values of the time t^ # 



Table V-V 

Comparison of Analytic and Computer Values of the Time t^ 
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Plant 3) 

The computer was set up as shown in Figure V-30 using the basic 
circuits of Section IV for 



(V-21) x 
Initial states of x(0) 



0 

-1 



to} . ra • - ® 



4 l 

were selected to 



provide systems requiring differing numbers of control signal switching. 

The criterion used for the selection of AT was A V) . • A Y ) . 0 

' l+l ' i 

for all three initial states. The data are presented in Table V-VI, 
the plots of >1 versus Y| are presented in Figures V-31, V-3 2, 
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and V-33 respectively for the three given initial states and the plots 
of X 2 (t) versus x-|_(t) are presented in Figures V-3U, V-3£, and V-36. In 
all three cases the convergence to v^° was rapid and x(t^) and 
(z(t^, rj ) + Xq) were of the same order of magnitude as 



Plant k) 

The computer was set up as shown in Figure V-37 using the basic 
circuits of Section IV modified for the three-dimension system 



(V-22 ) x 



0 10 
0 0 1 
0 0 0 



u , |u < 1 



ix(0) = ) 1 r was the initial state which was investigated. Initially 
the criterion for the selection of AV was that A"v^ + ^ • A^J ^ 
however, after seven iterations it was decided to use as the criterion, 



AtT 



which would yield t- 



li+1 



> t l. 



Using this criterion, a total 



of sixty iterations were completed at which time x(t^) 0 and 

t(t ± ) -*(0). The data for iterations fifty through sixty are 

presented in Table V-VII. There was no apparent pattern for the selection 
of AT'. In all iterations it was possible to duplicate values of t^ 
for successive re-iterations with the same comouter parameters; it was 
not possible to duplicate values of zT(t^). With this in mind some error 
analysis work is indicated to determine the effect of an error in z(t^), 
introduced by some random signal, has on the convergence to ^ 
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VI 



SUMMARY OF RESEARCH RESULTS 



The optimum control law for two-dijnension systems is readily 
obtained using Neustadt's synthesis technique as adapted to analog 
computer solution* For the solution of three-dimension and higher 
order systems some refinement of the computer technique used in this 
research is necessary. 

The initial choice of >7 has relatively little effect on the 
convergence to the optimum steering order. There are a wide variety 
of criteria for the selection of ATr which in general lead to 
satisfactory results for the two-dimension system. For design 
of an automatic system perhaps the criterion that A • A rj ^ > 0 

would be the best suited. In general large values of ATT give slow 
convergence with large oscillations about v}° while small values of 
ATT give slow positive convergence. There obviously are some At" 
which give the best convergence but no optimization was attempted 
in this research work. 

Using the finite difference equation in computing iterative 
values of Yj it becomes possible for vj to leave the domain. 
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VII. POSSIBLE EXTENSIONS OF THIS RESEARCH WORK 



1. The natural extension of this work is to continue on into 
the third and higher dimension systems. 

2. A significant extension of this work would be the optimization 
of AT. 

3. An error analysis of the effect random errors in the computation 
of z(t^) has on the convergence to yf 0 would be noteworthy. 

i|. To be of maximum use the system must be fully automatic, 
utilizing logic circuits, and/or a digital-analog combination which 
subject requires some investigation. 
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